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Abstract. It has been known for a long time (Mes- 
tel 1953) that the meridional circulation velocity in stars, 
in the presence of ^-gradients, is the sum of two terms, one 
due to the classical thermal imbalance (fi-currents) and 
the other one due to the induced horizontal ^-gradients {fi- 
induced currents, or /i-currents in short). In the most gen- 
eral cases, /^-currents are opposite to fi-currents. Simple 
expressions for these currents are derived under some sim- 
plifying physical assumptions presented in the text, and 
their physical interpretations are discussed. Computations 
of the and //-currents in a 0.8 M0 halo stellar model in- 
cluding classical element settHng show that the /i-currents 
are larger (in absolute values) than the fi-currents in all 
the star: some new physics has thus to be invoked in this 
case. 

We show here how such processes could possibly lead 
to a quasi-equilibrium stage in which both the circulation 
and the helium settling could be cancelled out. As lithium 
diffuses in the same way as helium, we expect a very 
small lithium concentration gradient below the convective 
zone in "plateau stars" (main-sequence Pop II stars) , much 
smaller than the one expected for pure element settling. 
This could possibly account for the very small dispersion 
observed for the lithium abundances at the surface of these 
stars. This should also have important consequences in 
other contexts which will be discussed in forthcoming pa- 
pers. 

The present computations show that element settling 
in slowly rotating stars leads to surface abundances which 
depend on the competition between //.-currents and fl- 
currents, in a way which had not been taken into ac- 
count in previous computations. This may change our gen- 
eral understanding of the diffusion processes of chemical 
species in rotating stars. 
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1. Introduction 

The present study was motivated by the so-called "lithium 
plateau paradox" (Vauclair 1999), which may be formu- 
lated in the following way: 

- From spectroscopic observations, the lithium abun- 
dance in main sequence Pop II field stars with effec- 
tive temperatures larger than 5500 K is remarkably 
constant. Since the first observations of the lithium 
plateau by Spite and Spite (1982), many abundance 
determination have confirmed the constancy of the 
lithium abundance in most of these stars, although a 
small slope may exist as a function of effective temper- 
ature (Thorburn 1994, Ryan et al. 1996, Spite et al. 
1996, Bonifacio and Molaro 1997). Moreover the dis- 
persion around the average value is extremely small. 
Whether this dispersion is real or below the observa- 
tional errors is still a subject of discussion. (Deliyannis 
et al 1993, Bonifacio and Molaro 1997, Molaro 1999). 
In globular clusters however, some dispersion exists for 
subgiants close to the turn-off (DeHyannis et al. 1995, 
Boesgaard et al. 1998). 

- On the other hand, theoretical computations lead to 
predictions of large variations of the lithium abun- 
dance from star to star, unless some "ad hoc" pro- 
cess forces the lithium abundance to remain con- 
stant (Charbonnel and Vauclair 1998, and references 
therein). The reason for these predicted variations is 
related to the diffusion processes which take place in 
the radiative regions inside the stars, below the outer 
convective zones. Due to pressure and temperature 
gradients, lithium settles down towards the stellar cen- 
ter, as well as helium and other heavier elements. Mix- 
ing induced by rotation, internal waves, or mass loss 
related motions may slow down the settling process, 
but then it brings up to the convective zone matter in 
which lithium has been destroyed by nuclear reactions. 
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These Htiiium depleting processes have been exten- 
sively studied in the litterature: (e.g. Michaud et al. 1984, 
Vauclair 1988, Pinsonneault et al. 1992, Vauclair and 
Cliarbonnol 1995 and 1998). The results are well known 
and can be summarized as follows: 

- If the radiative zones were completely stable, lithium 
would be depleted in Pop II stars due to settling, 
and the effect would increase with effective tem- 
perature. As a consequence, the lithium abundances 
should decrease with increasing Tes in contradiction 
with the observations. This result suggests that some 
macroscopic motion must occur and compete with the 
lithium settling, (e.g. Michaud et al.l984, Vauclair 
1988, Chaboyer et al. 1992). 

- Rotation-induced mixing, as prescribed by Zahn 
(1992), or Pinsonneault et al. (1992), could account 
for the observations. However the small dispersion of 
the observed abundances is difficult to interpret in this 
framework (Vauclair and Charbonnel 95) and gives in 
any case an upper limit on the initial lithium value 
(Pinsonneault et al. 1998). 

- Mass-loss could also be a possible explanation (Vau- 
clair and Charbonnel 1995). However all the stars 
should have suffer exactly the same mass loss rate, 
about 30 times the solar wind rate, during all their 
Hves. 

In all cases, the lithium plateau can be accounted for 
only with "ad hoc" hypothesis, where some parameter is 
assumed constant (or with very small variations) in all 
stars (mass loss rate, rotation rate...). This is not satisfy- 
ing, and observers often prefer to forget about this physics 
and decide that the lithium abundance observed in Pop II 
stars must be the primordial one. 

This interpretation however is in contradiction with 
the strong improvements of stellar physics obtained these 
last few years. The theory of stellar structure and evolu- 
tion includes the best available plasma physics, with new 
equations of states, opacities, nuclear reaction rates and 
element settling. Helioseismology represents an excellent 
tool for testing this physics: the agreement of the sound 
velocity in the models and in the "seismic Sun" (deduced 
from helioseismic modes) is much better when element set- 
tling is introduced (e.g. Bahcall and Pinsonneault 1995, 
Richard et al. 1996, Brun et al. 1998, Vauclair 1998). 

This is the reason why the lithium plateau becomes 
a paradox, with observations in contradiction with the 
predictions of the best available theory. What does the 
lithium plateau want to tell us that we have not yet un- 
derstood ? 

The clue may be related to the physics of rotation- 
induced mixing in the presence of /x-gradients. This has 
been studied by several authors, beginning with Mestel 
(1953, 1957, 1961, 1965) and Kippenhahn (1958, 1963), 
and recently by Zahn (1992) and Maeder and Zahn (1998). 



In the present paper, the physics of meridional circulation 
in the presence of /x-gradients is revisited. We show that, 
when one takes into account helium gravitational settling, 
/x-induced currents opposing the traditional meridional 
currents (hereafter called fi-currents) take place, and that 
they can become of the same order of magnitude. Then 
the circulation freezes out. Furthermore, in this case the 
element settling itself may also be prevented. We are in 
the presence of a self-regulating process where the rota- 
tional currents and the hehum settling cancel each other. 
As lithium diffuses in a way similar to helium (at the same 
rate), the lithium settHng is also cancelled by this process. 
We show how this could account for the small dispersion of 
the lithium plateau. It may also have other consequences 
which will be studied in forthcoming papers. 

In section 2, the theory of meridional circulation in- 
cluding the effects of /x-gradients is recalled. Orders of 
magnitude for the case of Pop II plateau stars are dis- 
cussed in section 3 and evidences of the importance of 
/x-induced currents are given. We show how element set- 
tling can be prevented when /x-induced currents cancel 
f2-currents. A general discussion is presented in section 4. 



2. Meridional circulation including /x-induced 
currents 

2.1. Classical meridional circulation 

In rotating stars, the equipotentials of "effective grav- 
ity" (including the centrifugal acceleration) have eUipsoi'- 
dal shapes while the energy transport still occurs in a 
spherically symctrical way. The resulting thermal imbal- 
ance must be compensated by macroscopic motions: the 
so-called "meridional circulation" (Von Zeipel 1924). The 
stellar regions outside the convective zones cannot be in 
complete radiative equihbrium. They are subject to en- 
tropy variations given by (see the review by Zahn 1993): 



pT 



dS_ 
'di 



. u • V = S ) = -V • F p£„ 

= pen (^ 0) 



(1) 



where F represents the heat flux, e„ the nuclear energy 
production and eo an energy generation rate which results 
from sources and sinks of energy along the equipotentials. 

The vertical component of the meridional velocitjf Ur 
is computed as a function of eq, in the stationary regime 
(from eq. 1): 



Ur = 



CppTj g 



which, for a perfect gas, reduces to: 



Vad 



9 Vad-V + V^ 



(2) 



(3) 
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where g represents the local gravity, Vad and V the usual 
adiabatic and real ratios ( — - — | and Vu the mean 



dlnP 

din II 



molecular weight contribution , ^ 

\ a m i-^ ^ 

The expression of en is computed by expanding the 
right-hand-side of eq. (1) on a level surface and writing 
that its mean value vanishes. 

For nearly-uniform rotation and negligible /^-gradients, 
one finds (cf Mestel 1965, Zahn 1993 and references 
therein): 



= -r TT 



3 M V GM 



2ttGp 



P2(cos6') (4) 



all quantities being computed at radius r. 

In this case, neglecting in (3), the vertical merid- 
ional velocity may be written: 



_ 8 _L Vad ^ 

" 3 M Vad-V 
Introducing Ur as: 



2-kGpJ 



Ur = Ur P2 (cos 6) 



P2(cos6») (5) 



(6) 



the horizontal meridional velocity is given by (see Mestel 
1965): 

U0 = - -J— (pr'^Ur)smdcosd (7) 
2pr dr 

In the presence of mean molecular weight gradients, 
these expressions must be modified, as discussed below. 



2.2. jjL-induced currents 

Mestel (1953) and (1957) (see also Mestel 1965) pointed 
out that, in the presence of vertical /x-gradients, eo con- 
tains other terms related to the resulting horizontal varia- 
tions of /i: the so-called "/i-induced currents". The expres- 
sion of Eq, obtained in this case, with the assumption of 
perfect gas law, has been derived in detail. 

More recently Maeder and Zahn (1998), hereafter 
MZ98, gave a modified expression for eq, which takes into 
account several effects which were not included in the pre- 
vious computations: more general equations of state in- 
stead of perfect gas law, presence of a thermal flux induced 
by horizontal turbulence, non-stationary cases. 

Here we wish to focus on the importance of non- 
negligible /x-gradients for the meridional circulation, which 
may modify our understanding of the low- rotation regimes 
in the presence of microscopic diffusion. 

The basic idea is that on a level surface (with constant 
pressure) specific quantities like temperature T, density 
p, gravity g, thermal conductivity x, rotation rate VI and 
mean molecular weight p, may vary around an average 



value. Following Zahn 1992 and MZ98 these quantities x 
are expanded on isobars as: 



x{P, 0) = x{P) + x{P)P2 (cos ( 



(8) 



The expansion of the energy term in eq. (1) leads to dcriva- 
d 

fives in — — which have to be transformed through an 

d 

equation of state to lead to the radius derivatives — . 

dr 

In this procedure, two kinds of terms finally appear: 

1) those which are directly related to the rotation rate, 
for example the gravity fluctuations 

2) those which are related to the local /x-gradient. 

As we discuss a basic physical point, we prefer to sim- 
plify the expressions under specific assumptions, neglect- 
ing for the moment the secondary terms. We assume per- 
fect gas law, and nearly rigid rotation (the importance 
of this assumption will be discussed in section 4). In this 
case, the gravity fluctuations are given by : 



9_ 
9 



GM ) 



(9) 



As we focus on the physical processes which occur be- 
low the convection zones in cool stars, we also neglect all 
the terms related to energy production (stellar cores will 
be discussed in a forthcoming paper). Keeping only the 
non-neghgible terms in MZ98's expressions, we obtain: 



en 



{Ea + E^)P2{cose) 



(10) 



with: 
En = 



0V\ 
GM ) 
r d 
3 dr 



2ttGp) 
Ht'^ - (x^ + XT + l)A 



(11) 



Here p represents the density average on the level surface 
(~ p) while pm is the mean density inside the sphere of 
radius r; Ht is the temperature scale height; A represents 

the horizontal a fluctuations ^; Xu ^^^d Yt represent the 

derivatives: 



dlnp, 



XT 



P,T 



91nx 
dlnT 



(13) 



p,^l 



In the following we will also neglect the "Gratton-Opik" 
term , which is justified below the convective zones 

in slowly-rotating cool stars. 



2.3. Physical interpretation 

In the presence of p-gradients, the circulation velocity is 
the sum of two terms: Eq, which we shall call the "Ct- 
current", and E^, , which we shall call the "/U- induced cur- 
rent" or, for short, the "/u-current". In the general case. 
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-E,, may also be related to fl through the horizontal /i- 
fluctuations A. We will see however that, in the stationary 
case, A is independent of O. In most of the star, Eq is pos- 
itive (except in the very outer layers) while E^^ is negative 
(except in case of strong second derivatives of A, which 
are supposed negligible in the following): the /i-currents 
are opposite to the fi-currents. When A varies directly 
like r (this situation occurs in the stationary case, as will 
be seen below), becomes: 



E„ 



-2 ^ 



Ht 



-A 



1 + 



(1 + XT + X^) 



r 



Which we will write: 



En — — 
3 



GM 



2pm Ht . 

~ -CA A 

p r 



(14) 



(15) 



(16) 



where ca is of order unity. 

The physical interpretation of these currents can be 
given in the following hand-waving way. Due to centrifu- 
gal effects, the effective gravity is constant on ellipsoidal 
shells (Figure 1) while the temperature gradients keep the 
spherical symetry. Consequently the temperature is larger 
at the poles than at the equator on a level surface, lead- 
ing to the f2-currents. The element abundances are also 
expected to varjf along level surfaces, as they depend on 
temperature (in case of gravitational settling, they also 
depend on gravity, so that they should not be constant on 
spheroids either). 

The density fluctuations on level surfaces are given by 
(Zahn 1992): 

1 r2 dVL^ 
3 g dr 



P 
P 



(17) 



As a consequence, in the approximation of nearly rigid 
rotation, level surfaces are both surfaces of equal pressure 
and density. In this case p, varies like T. 

In first approximation, all other parameters assumed 
constant, the specific entropy is proportional to the num- 
ber of particules per unit volume (Sackur- Tetrode formula, 
see Vauclair 1993), or inversely proportional to p. Thus 
the horizontal /i-gradients lead to a larger entropy at the 
equator than at the pole, which explains why /x-currents 
are opposite to fi-currents. 

The /Lt-currents induced by the flattening of the level 
surfaces are very small though, so that they do not prevent 
the onset of the meridional circulation. However, as soon 
as the circulation begins, the horizontal /x-gradients in- 
crease as the upward flows bring /x-enriched matter at the 
poles while the downward flows bring /i-dcplctcd matter 
at the equator. As will be seen below, a stationary solu- 
tion can occur for which, in some cases, the /U-currents 
may become of the same order as the 0-currents, leading 
to new interesting physics. 




Fig. 1. Schematic drawing of a "level surface" in a rotating 
star. Temperature and mean molecular weight are larger 
at the poles than at the equator, inducing a competition 
between two effects: the so-called "f2-currents" from pole 
to equator and "^-currents" from equator to pole. At equi- 
librium, the circulation stops. Then element settling pro- 
ceeds at the pole in a more rapid way than at the equator, 
thereby decreasing the horizontal /i-gradient and Ep. The 
circulation begins again with matter going up at the pole: 
it brings back the material in which helium has settled, 
restablishing the equilibrium horizontal /i-gradient. Such 
a process can cancel altogether meridional circulation and 
helium settling. 

2.4- Discussion on the horizontal variations of p 

The order of magnitude of the horizontal /i-gradients in- 
duced by the flattening of the level surfaces A/ can be 
obtained with the approximation that the radius differ- 
ences between the polar and equatorial regions are given 
by: 



A = 



Ar _ fiV^ 
~V ~ GM 



(18) 



The values of this flattening parameter are given in Fig- 
ure 2 as a function of r in a 0.8 Mq halo star which rotates 
with a velocity of 10 km.s~^. 

In this case, the horizontal p fluctuations can be writ- 
ten: 



p 

A/ = - ~ • Vln/z 



(19) 



p GM 

or, with the following deflnition of the "^-scale height": 

H^ = {y\np)-^ (20) 

the horizontal fluctuations related to the flattening be- 
come: 



A, 



2^3 



GM H„ 



(21) 
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These fluctuations are small compared to those which oc- 
cur as soon as the meridional circulation is established. In 
case of a laminar circulation (pure advection) the order of 
magnitude of A is given by: 



Aa ~ r ■ V In /I ~ — 



(22) 



According to Zahn (1992) and Chaboyer and Zahn (1992), 
the shears induced by the circulation lead to a strong 
horizontal turbulence, which reduces the horizontal /i- 
gradients. In this case, introducing a horizontal diffusion 
coefficient Dh , the local variations of /i are solutions of 



9/i 
'dt 



dr 



(23) 



In the stationary case, we obtain for the horizontal /i- 
gradient including turbulence: 



A 



_ ^ = 



dr 



6D. 



(24) 



With the assumption that Dh is proportional to Ur-r and 
not simply to Ur as written by mistake in MZ98, eq. (22) 
becomes: 



A. 



1 9 In /i 
6ah 9 In r 



1 r 



Replacing in eq. (16) gives the simple expression: 



^ 3a/i i/^ p 



(26) 



3M 



with p ^ p and pm = 

In summary, in the most likely case of meridional cir- 
culation associated with a strong horizontal turbulence, 
the vertical velocity of meridional circulation can simply 
be written: 



with Ur 



En - 



f/^ Pa (cos 6*) 

1 / Vad 

.9 " 



Vari - V 



3GM 
Ca M Ht 





^ ah 4:TTr^p Hj 



(28) 
(29) 



According to Zahn's (1993) prescriptions, while the hor- 
izontal diffusion coefficient is larger than Ur ■ r {Dh = 
ahUr ■ r with ah > 1), the vertical transport may be ap- 
proximated as an effective diffusion coefficient of order: 



D. 



'off 



C, 



-Ur-r with Ch — 30 



^ Following Chaboyer and Zahn (1992) and MZ 98, is de- 
velopped on r instead of P. The difference is negligible though, 
as the flattening effect on fi is extremely small. 




flattening parameter 
V(rot)= lOkm.s'' 



0.4 0.6 

r/R 

Ar il'^r^ 

Fig. 2. Flattening parameter A = — ~ "gm" '^^ ^ ^'^ ^® 
plateau star with Vrot — 10 km.s^^ (see text for details). 
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Fig. 3. Helium abundance profile and grad(ln/i) in a 
0.8 M0 plateau star, including element settling, at the age 
of 12 Gyr. The increase in grad(ln/i) below the convective 
zone is clearly related to the "^He slope. 



The reason why i?^ does not depend on Q but depends 
on may be physically interpreted in the following way. 
Suppose that at level r the meridional circulation can be 
approximated by an upwards and a downwards moving 
fiows. The horizontal turbulence, described by the hori- 
zontal turbulent diffusion coefficient Dh, leads to matter 
transfer from one flow to the other. In this framework, a 
vertical /.t-gradient induces a horizontal /^-gradient which 
plays a non-negligible role in the entropy balance and has 
a feed-back effect on the velocity Ur- Increasing the rota- 
tion velocity should increase Ur, the horizontal //-gradient 
and the corresponding Dh. But the increase of the mixing 
effect related to Dh decreases this horizontal /(-gradient. 
The horizontal turbulence "swallows" the possible varia- 
tions induced by different rotation velocities, leaving the 
vertical mixing process unchanged. 
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Fig. 4. Comparison of the /i-gradient 



with the 



dlnP 

usual adiabatic and radiative gradient differences. The /Li- 
gradient is always negligible compared to the others, ex- 
cept in the core and right at the bottom of the convective 
zone. 

3. Orders of magnitude and discussion 

3.1. Comparison between the ^.-currents and the 
fjL-currents 

At this point, it is interesting to compare the orders of 
magnitudes of Ei^i and in stellar models. As the present 
paper has been motivated by the existence of the lithium 
plateau in halo stars, we chose to discuss as an example 
the case of a 0.8 M© low metallicity star, situated at the 
hot end of the plateau {Z = 10~^; age: 12 Gyr). Figure 2 

gives the values of the flattening parameter A(r) = , 
with the assumption of a stellar rotation velocity Kot = 
10 km.s"^ which represents an upper limit for these slowly 
rotating stars. Figure 3 displays the abundance profile in 
the same star, including clement settling (Charbonncl and 
Vauclair 1999). It also shows the yU-gradient Vhiji, which 
is essentially related to the helium profile. 

The values of Vu = ( I are given in Figure 4 

\omF J 

also for the model including element settling, together 
with the difference Vad — V: we can see that is al- 
ways negligible except in the central core and right at 
the bottom of the convective zone. The /x-scale heights 



H^, defined as 



dr 



are given in Figure 5 for two 



models of the same star, one without element settling (st) 
and the one including element settling (dif ) . The tempera- 
ture scale height Ht, computed in the model with element 
settling, is also shown (the difference in Ht between the 
two models is negligible) . The effect of helium diffusion is 
clearly visible on this graph. 

Figure 6 presents for these two models a comparison of 
the absolute values of -B^ and Eq. Here again the rotation 
velocity is taken as 10 km.s~^ which represents an upper 
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Fig. 5. /i-scale heights computed in two different 0.8 Mq 
plateau stars models: one with no element settling in- 
cluded (st) and one with element settling included (dif). 
The temperature scale-height is also shown for the model 
with element settling (there is no significant difference 
with the other model). The difference in the /U-scale 
heights is the signature of helium settling. 

limit for halo stars. The Efj, values are computed with 
the assumption of meridional circulation and horizontal 
turbulence, and the diffusion coefficient Dh is assumed 
equal to 10 Ur-r {an = 10). We can see that in the model 
without settling, is larger than in the central 
regions and smaller in the outer layers while in the model 
including settling, \E^\ is always larger than \Efi\. This 
means that some physical process must occur which has 
not yet been included in the theory. The case of the stellar 
core, where other effects can occur, will be discussed in a 
forthcoming paper. Here we focus on the radiative region 
below the convective zone. 

Suppose the star begins on the main sequence with 
homogeneous abundances. Then \E,,,\ lies below \Eu\ and 
the meridional circulation can occur in the same way as 
we have already discussed. Horizontal turbulence may de- 
velop, and a vertical mixing takes place with an effective 
diffusion coefficient Deff (Zahn 1992). This vertical mix- 
ing slows down the element settling without stopping it. 
Defining He as the concentration scale height, the element 
settling would stop only for: 



Hr_ 



H„ 



where Dm is the microscopic diffusion for the considered 
element and k^, the coefficient in front of the pressure gra- 
dient in the microscopic diffusion velocity (see, for exam- 
ple, Vauclair and Vauclair 1982). However, before reaching 
this equilibrium concentration, the induced /i-currents 
become of the same order as the f2-currents, there- 
fore stopping the circulation. 

We can compute, for the same 0.8 Mq star, the critical 
/Lt-scale height for which \Ei^\ = \Eq\. 
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Fig. 6. /i-currents (in absolute values) computed in the 
two 0.8 M0 plateau stars models, with and without ele- 
ment settling, and fi-current computed in the model with 
setthng, assuming Vrot = 10 km.s~^ and a — 10 (see text). 
We can see that in the model with setthng, is always 
larger than \Eq\. As most halo stars have rotation veloci- 
ties below this value, \Eq\ is expected to be even smaller 
than shown here, which reimforces this conclusion. 



We obtain: 



Ht Pm 

8a/,A p 



ah 32irn'^pr^ 



(31) 



The profile of H™'^ (r) is given in Figure 7, also for T^-ot = 
10 km.s"^ and ah ~ 10. For different rotation velocities 
H"'^ is simply multiphed by (100/wrot)^. 

It is interesting here to compare this critical /U-scale 
height with the concentration scale height given by equa- 
tion (30) . For a simple mixture of completely ionised hy- 
drogen and helium, the mean molecular weight is given 
by: 



2 + 3c 



With c of the order of 0.1, we find for helium: 



(32) 



(33) 



Equation (30) gives, for the equilibrium concentration 
scale height: 



Hr = H, 



Taking kp ^h,Hp~ 10^° and D^s/Dm ^ 100 we find: 



-Peff 

'^k D 



(34) 



HrC^2x W cm and H„ 



10^^ cm. 



We can check from Figure 7 that the critical H^j which 
stops the circulation is reached much before the equilib- 
rium He which would stop the helium settling. 

Now the question arises: what happens when be- 
comes equal to \Eq\ ? 
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Fig. 7. Critical /K-scale height for which \Efj\ = \Eq\ in 
a 0.8 M© plateau star, with the same assumptions as in 
Figure 6. In this case, the /i-scale height is more than 5 
orders of magnitude larger than the stellar radius: the re- 
sulting helium (and also lithium) profiles should be fiat be- 
low the convection zone. For lithium, it should remain fiat 
down to the nuclear destruction region. As -ff"'*, varies 
like (tirot)"^, this conclusion is reimforced for slower rota- 
tions. 



3.2. Self -regulating process in case of meridional 
circulation and settling 

Let us first summarize the situation of a slowly rotating 
star in which element settling leads to an increase of the 
/i-gradient below the outer convection zone. 

At the beginning, the star is homogeneous and merid- 
ional circulation can occur, leading to upward fiows in the 
polar regions and downward fiows in the equatorial parts 
(except in the very outer layers where the Gratton-Opik 
term becomes important, which we do not discuss here). 
The /x-currents, opposite to the classical fl-currents, are 
first negligible. The /i-gradients increasing with time be- 
cause of hehum setthng, the order of magnitude of the 
/i-currents also increase and the /U-scale height decreases 
until it reaches H'r"*' for which the circulation vanishes. 

This does not occur all at once: as H"^"^ decreases with 
depth below the convective zone (Figure 7) we expect that 
the meridional circulation freezes out step by step. An 
equilibrium situation may be reached, in which the tem- 
perature and mean molecular weight gradients along the 
level surfaces are such that f2-currents and /U-currents can- 
cel each other. 

Once it is reached, this equihbrium situation is quite 
robust. Suppose that some mechanism leads to a decrease 
of the horizontal //-gradient: then \E^\ becomes smaller 
then \Eii\ and the circulation tends to be restablished in 
the \Eq\ direction. The fiow goes up in the polar regions, 
bringing back matter with larger p.: the p gradient is re- 
stored and the circulation stop. 
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Suppose now that the horizontal ^-gradient is in- 
creased. Then jE"^! becomes larger than |£^n| and the cir- 
culation begins in the direction. Matter with larger ji 
goes up in the equatorial region, decreasing the /i-gradient. 
Here again the circulation stops. 

As there is no more circulation, we would expect that 
helium settling can proceed further. The microscopic dif- 
fusion velocity, in case of pure gravitational settling, can 
simply be written (Vauclair and Vauclair 1982): 

Vd = -Dra^ = -Dmkp^ (35) 

with: 

D a T^/^/p (36) 

On a level surface, fi varies like T, so that Vd a T^/^ and 

AvD _ 5 AT _ 5^ 

VD ~ 2 ~ ~ 2 

As a consequence, on a level surface, the helium gravi- 
tational settling is more efficient in the regions of larger T, 
namely the polar regions, than in the equatorial regions. 

We can now describe the scenario: once the merid- 
ional circulation is frozen, the helium gravitational settling 
would like to proceed and lead to a decrease of the hori- 
zontal /i-gradient. However the /i-currcnts become smaller 
than the fi-currents and /i-enriclied matter is brought back 
into the polar region, restoring the original /i-gradient. 

We may be here in the presence of a self-regulating pro- 
cess in which the /x-gradicnts which should increase with 
time due to the settling are prevented to do so because of 
the currents equilibrium. Such a process would stop alto- 
gether the meridional circulation and the element settling 
! 

In this case the clement abundances that we observe 
at the surface would be the result of the self-regulating 
process and not the result of pure settling. For halo stars, 
as shown on Figure 7, this would occur for iJ^ much larger 
than the stellar radius: the abundance profiles would be 
flat below the convection zone. 

This process could be the reason for the low dispersion 
of the lithium abundance in the lithium plateau of halo 
stars: the lithium abundance profile would be fiat down 
to the place where the diffusion time scale is of the or- 
der of the nuclear destruction time scale, that is the place 
of the lithium peak (Vauclair and Charbonnel 1998). The 
observed abundance would then depend only on the mi- 
croscopic physical parameters, constant from star to star, 
which would explain the very small observed dispersion. 
This process will be studied numerically in a forthcoming 
paper. 

4. Conclusion 

Computations of meridional currents in the presence of 
/i- gradients show that, as already pointed out long ago by 
Mestel (1953), and recently studied by Zahn (1992) and 



Maeder and Zahn (1998), //--induced currents (/^-currents 
in short) settle in opposition to fi-currents. In the present 
paper, we have computed the orders of magnitude of these 
currents under simplifying assumptions. Among these as- 
sumptions (energy production negligible, perfect gas equa- 
tion of state, Gratton-Opik term negligible, nearly solid 
rotation), the most stringent one is the assumption of 
nearly-solid rotation, as pointed out by Zahn (1999). As- 
suming solid rotation means that some other process like 
internal waves must have forced the transport of angular 
momentum (Zahn et al. 1997). In any case, we know from 
helioseismology that the Sun does indeed rotate as a solid 
body below the outer convcctive zone (Brown et al. 1989). 
This is not yet quite understood, but we may infer that 
the same process which forces solid rotation inside the Sun 
also acts in other solar-type stars. This may not be crucial 
for our general conclusions, although including differential 
rotation would need a more complicate treatment of the 
equations, and would add another unknown parameter. 
This may be studied in the future. 

We have shown that under these assumptions /U- 
currents may become of the same order of magnitude as fl- 
currents in case of element settling, long before the concen- 
tration gradient induced in case of pure settling is reached. 
In the case of a 0.8 M© halo star with standard helium 
settling, we have computed that the /x-currents would be 
several orders of magnitude larger than the 0-currents ! 
This is also the case in the central stellar regions, due to 
nuclear reactions. 

Here we have only focused on the physical processes 
which occur below the outer convective zones. In the case 
of halo stars we found that, for slow rotation, /^-currents 
cancel 0-currents for very small concentration gradients, 
of order 10''^^ (Hf^ = 10^^ cm). Furthermore, in this 
case helium setlling may also be cancelled out as it would 
decrease the effect of horizontal /i-gradients, thereby set- 
ting again meridional currents which would restore the 
original gradients. This may become a self-regulating pro- 
cess which could explain why the observed lithium in the 
plateau stars has such a small dispersion: the inferred 
lithium gradient inside the star, below the convective zone, 
would be very small and the lithium profile nearly fiat un- 
til the lithium nuclear destruction region. This would give 
a surface lithium abundance close to that of the "lithium 
peak" (Vauclair and Charbonnel 1998), the same one for 
all the stars. 

Numerical simulations of the whole process have to 
be done, to check all these effects. It will also be tested 
for the solar case (including central regions), and for the 
case of Am stars (for which the Gratton-Opik term is not 
negligible), in forthcoming papers. There are many obser- 
vations in stars which give evidences of mixing processes 
occuring below the outer convection zones as, for example, 
the lithium depletion observed in the Sun and in galactic 
clusters. The process we have described here should not 
apply in all these stars. The reason could be related to 
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the rapid rotation of young stars on the ZAMS and to 
their subsequent rotational braking (here we suppose that 
the stars always rotated slowly). Also for stars in which 
the Gratton-Opik term is no more negligible, the whole 
process has to be revised. Here we have only studied equi- 
librium situations, while time-evolving processes should 
be tested in the future. Finally other mixing processes, 
like ABCD or GSF instabilities, could play a role. 

In any case the present computations show that el- 
ement settHng in slowly rotating stars leads to surface 
abundances which depend on the competition between ^- 
currents and fi-currents. This may change our general un- 
derstanding of the diffusion processes of chemical species 
in rotating stars. 
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